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SUMMARY 


We  present  computational  test  results  from  a  double 
precision  version  of  the  FORTRAN  program  'WAGNER'.  The 
potential  value  of  this  integral  equation  technique  is  its 
purported  ability  to  model  propagation  losses  over 
two-dimensional  paths  through  irregular  and  inhomogeneous 
terrain..  However  the  numerical  integration  algorithm  can 
be  affected  by  poor  convergence  and  inaccuracy, 
partr.culcrly  at  frequencies  above  the  MF  band.  We  have 
investigated  these  properties  in  terms  of  selected  key 
parameters  and  have  derived  empirical  criteria  governing 
convergence  and  physical  accuracy.  We  find  that 
calculations  of  basic  transmission  loss  using  WAGNER  are 
subject  to  significant  limitations  at  HF  frequencies  as 
the  terrain  becomes  irregular  and  as  the  conductivity  and 
dielectric  constant  of  the  ground  both  become  smaller.- 
On  the  other  hand,  the  shadow  loss  due  to  an  obstacle  may 
be  reliably  calculated  by  WAGNER  on  the  assumption  that 
shadow  loss  is  independent  of  ground  electrical  constants, 
thus  permitting  the  use  of  dummy  values  of  conductivity 
and  dielectric  constant  to  ensure  that  the  physical 
accuracy  criterion  is  satisfied. 
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1 .  INTRODUCTION 

Ott(re£.l,2)  has  formulated  an  integral  equation  technique  for  the  calculation 
of  ground  wave  propagation  losses  over  paths  characterised  by  irregular 
terrain  and  varying  electrical  ground  constants  (conductivity  and  dielectric 
constant).  The  calculation  is  performed  by  step-wise  numerical  integration 
along  the  transmission  path  of  the  sum  of  a  c  irect  ray  and  a  ground-reflected 
ray,  and  is  coded  in  the  FORTRAN  program  'WAGNER ’ . 

References  1  and  2  presented  example  calculations  at  frequencies  between  1  and 
10  MHz  showing  results  which  are  consistent  with  those  obtained  using 
alternative  methods.  Generally  favourable  comparisons  with  selected 
experimental  measurements  were  also  presented  in  references  3  and  4, 
indicating  the  reliability  of  the  technique.  The  technique  is  however  subject 
to  certain  qualifications,  for  example: 

(a)  the  technique  is  spatially  two-dimensional  and  treats  only  propagation 
along  the  great  circle  path  betwean  two  antennas,  neglecting  side-  and 
back-scattered  radiation;' 

(b)  csre  must  be  taken  to  ensure  that  the  integration  step  length  is 
sufficiently  small  to  enable  the  solution  to  converge;  hovevei  too  small  a 
rtep  length  will  lead  to  excessively  large  pregram  execution  times; 

(c)  a  frequency-coupled  limitation  exists  on  th*  maximum  gradient  of  path 
terrain  which  may  be  treated  accurately;  tnis  "heuristic  uncertainty 
principle" (rel. 1)  is  given  as 


f  <  1^ 

<\J 


where  y1  is  the  gradient  and  f  is  the  frequency  in  MHz,  and  shows  that  a 
gradient  of  0.25  should  return  accurate  loss  calculations  throughout  the  HF 
band  up  to  a  maximut  frequency  of  4G  MHz.-, 

(d)  accuracy  can  be  machine  dependent,  with  small  wordlength  machines 
increasing  the  susceptibility  to  round-off  errors  unless  double-precision 
floating  point  arithmetic  is  used. 

(o)  computations  are  subject  to  poorly  understood  numerical  instability 
under  certain  conditions . 

However,  the  potential  advantages  of  the  technique  appear  to  outweigh  the 
disadvantages.  Within  the  constraints  of  the  '’heuristic  uncertainty 
principle"  it  should  be  possible  r.o  treat  a  path  having  an  arbitrary  terrain 
profile  and  include  variations  in  ground  electrical  constants,  for  example 
across  a  coastline.  Reliability  in  such  a  calculation  could  permit  a 
significant  reduction  in  the  need  to  perform  propagation  experiments  in  the 
field,  with  consequent  savings  of  time,  expense  and  manpower-  While  the 
residue  series  expansion  method(ref.5-6)  provides  a  reliable  and  quirk  means 
of  calculating  losses  over  smooth,  homogeneous  terrain  it  is  unable  to 
accommodate  terrain  irregularities,  in  contrast  to  WAGNER. 

Alternative  techniques  for  the  estimation  of  losses  due  to  terrain 
irregularities  are  based  on  knife-edgn  diffraction.  Bui lington(ref .7 < 
presents  a  chart-  from  which  it  is  possible  graphically  to  estimate  the 
attenuation  caused  by  c  single  knife-edge  located  between  two  antennas,  in  free 
space..  Its  adaptation  to  the  case  of  a  single  knife-edge  or.  the  surfac  •  of  a 
smooth  spherical  Earth  has  resulted  in  a  rhart  which  gives  the  shadow  less  of 
the  obstacle  with  respect  to  smooth,  spherical  homogeneous  Earth  ground  wave 
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propagation,  and  has  been  verified  experimental ly( ref .8) .  Multiple 
knife-edges  are  not  accommodated,  however,  and  must  be  treated  by  other 
techniques,  such  as  WAGNER. 

We  point  out  that,  while  s  hill  may  be  conveniently  modelled  either  by  a 
Gaussian  profile  or  by  a  knife-edge  for  the  purpose  of  calculating  wave 
attenuation,  results  differing  by  several  decibels  are  likely  to  be  obtained 
from  the  two  methods,  particularly  in  close  proximity  to  the  hill.  However, 
the  results  of  Ott(ref.l)  indicate  that  such  discrepancies  diminish  with 
distance  and  becoma  unimportant  at  distances  exceeding  on  the  order  of  5  times 
the  1/e  full  width  or  the  hill.  We  therefore  believe  that  with  this 
qualification,  comparisons  between  Gaussian  and  knife-edge  profiles  may  be 
made. 


1 . 1  The  program 

The  mathematical  details  concerning  the  integral  equation  for  ground  wave 
propagation  over  irregular  and  inhomogeneous  terrain  have  been  described  by 
0tt(ref.l,2)  and  will  not  be  repeated  here.  The  solution  to  the  resulting 
linear  Volterra  integral  equation  is  solved  by  the  method  of  Wagner(ref .9) , 
whence  the  title  of  the  computer  program.  This  solution  involves 
subdivision  of  the  transmission  path  into  steps  of  the  order  of  a 
wavelength  in  size  and  a  summation  in  signal  strength  over  all  subintervals 
along  the  path.  The  techniques  involved  in  obtaining  the  numerical 
solution  are  r.on-tiivial  and  at  this  time  a  formal  error  analysis  of  the 
procedure  apparently  has  not  been  made.  It  is  certainly  beyond  the  scope 
of  this  report.  Consequently  the  user  of  WAGNER  faces  a  lack  of 
information  on  such  matters  as  the  convergence  properties  of  the  technique 
and  the  accuracy  with  which  it  models  transmission  losses,  particularly 
within  the  HF  band. 

Our  version  of  WAGNER  was  originally  coded  in  single-precision  but  has  been 
converted  to  double-precision  to  reduce  round-off  errors..  Named  W5D,  it  is 
compiled  under  VAX  FORTRAN  version  4.5  and  runs  on  the  HF  Radar  Division 
VAX  8200  under  the  operating  system  VMS  version  4.4.  We  have  tested  the 
program  for  the  case  presented  in  figure  1  of  reference  1,  that  is,  a 

frequency  of  1  MHz,  ground  conductivity  of  0.01  S  m  1 ,  dielectric  constant 
of  10,  an  otherwise  flat  path  of  length  15  km  containing  a  Gaussian  hill  of 
height  1  Iud  located  at  a  distance  of  5  km  from  the  transmit  antenna.  The 
integration  subinterval  cr  step  size  used  in  reference  1  is  not  specified 
but  our  work  employed  s  step  size  of  100  m.  The  antennas  are  vertically 
polarised  and  located  a>,  ground  level.  The  resulting  basic  transmission 
loss  is  graphed  as  a  function  of  distance  along  the  path  in  figure  1(a). 
Agreement  with  reference  1  is  excellent,  validating  our  version  of  the 
program  at  this  frequency  in  the  MF  band.  Note  that  the  data  in 
reference  1  are  presented  in  terms  of  the  amplitude  attenuation 
function  (f)  with  respect  to  twice  the  free  space  field  which  is  related  to 
che  amplitude  basic  transmission  loss  at  wavelength  \,  between  isotropic 

antennas  separated  by  a  distance  x.  according  to 


r  — 
%  " 


2f 


X 


Our  calculations  are  presented  in  terms  of  basic  transmission  loss  in 
decibels  (dD).  LosSba  (»re  represented  by  negative  .lumbers.  The  execution 
tine  required  to  produce  the  results  of  figure  1(a)  was  250  s.  In  general, 
execution  time  varies  as  the  square  of  the  number  of  integration  steps. 
We  point  out  that  at  15  km  the  shadow  loss  resulting  from  the  hill  is  about 
-10  dB  with  respect  to  smooth,  sphecical  Earth  propagation,  in  excellent 
agreement  with  the  graphical  estimate  from  reference  8.  However  if  we  move 
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out  of  the  MF  band  into  the  HF  band  we  encounter  problems.  Figure  1(b) 
shows  results  for  a  frequency  of  10  MHz,  all  other  parameters  remaining 
unchanged.  It  indicates  that  the  shadow  loss  of  -10  dB  at  1  MHz  has  at 
10  MHz  been  transformed  into  a  gain  of  +41  dB.  Since  in  this  case  the 
knife-edge  technique  graphically  predicts  a  shadow  loss  of  -18  dB,  it  is 
clear  from  the  discrepancy  of  59  dB  that  WAGNER  has  broken  down. 


2.  TEST  TECHNIQUES 

Our  initial  attempts  at  applying  WAGNER  to  real  terrain  at  HF  frequencies 
produced  results  such  as  figure  1(b)  which  were  obviously  in  error.  In  view 
of  the  difficulties  involved  with  a  formal  error  analysis  it  was  decided  ‘.o 
carry  out  an  "experimental"  study  of  the  intervals  in  certain  key  parameters 
for  which  the  results  of  WAGNER  could  be  relied  upon.  The  aim  was  to  develop 
empirical  or  heuristic  "rules  of  thumb"  which  could  be  applied  to  predict 
whether  a  proposed  calculation  woulu  be  successful  or  not. 

The  key  parameters  selected  for  this  study  were  integration  step  size, 
frequency,  path  length,  terrain  slop.',  ground  conductivity  and  dielectric 
constant.  The  question  of  the  reliability  of  any  calculated  result  has  two 
aspects.  Firstly  it  is  necessary  to  know  whether  the  calculation  has  arrived 
at  the  convergence  limit  in  the  mathematical  sense.  Since  the  program  does 
not  iterate  on  the  integration  step  size  by  testing  the  difference  between 
successive  loss  calculations,  it  is  not  possible  to  know  "a  priori"  how  far 
any  given  result  deviates  from  its  convergence  limit.  Secondly,  successful 
mathematical  convergence  of  an  algorithm  does  not  necessarily  produce  a 
correct  answer  from  the  physical  viewpoint.  It  is  therefore  necessary  to 
compare  each  convergent  solution  with  the  loss  estimated  using  an  independent 
method  for  example  the  residue  series  expansion  algorithm,  in  the  case  of  a 
smooth  Earth,  or  the  graphical  knife-edge  method  in  the  case  of  a  single 
terrain  obstacle.  Thus  we  aim  to  produce  two  criteria,  one  describing  the 
conditions  necessary  for  convergence  and  the  other  similarly  specifying 
physical  accuracy. 

Both  convergence  and  accuracy  were  mapped  at  regular  grid  points  in  the  space 
of  the  integral  equation's  key  parameters.-  Calculations  of  basic  transmission 
loss  were  carried  out  with  uniformly  spaced  integration  steps  ranging  down  in 
size  from  1000  to  10  m,  at  frequencies  between  1  and  30  MHz,  path  lengths  from 
10  to  100  km,  terrain  gradients  from  0  to  3,  ground  conductivities  from  0.001 

to  0.1  S  m  1  and  dielectric  constants  from  1.1  to  100. 

Two  types  of  terrain  were  employed,  a  smooth  spherical  homogeneous  Earth  of 
effective  ridius  factor  equal  to  4/3  and  a  single  Gaussian  hill  based  on  the 
terrain  of  figure  1  but  permitting  variable  hill  heights  and  ground  electrical 
constants.  Note  that  the  maximum  gradient  of  terrain  is  given  in  this  case  by 

approximately  10  *  times  the  hill  height  in  metres. 


3 .  RESULTS 


3 . 1  Convergence 

3.1.1  Effect  of  path  length  and  frequency 

Basic  transmission  losses  were  calculated  over  a  sffioth  spherical  Earth 
at  path  lengths  L  =  10,  20,  50  and  100  km  using  integration  step  sizes 
Ax  ranging  from  1000  m  down  to  10  m.  A  ground  conductivity 

o  =  0.003  S  ra  1  and  dielectric  constant  K  =  3  were  used  at  each  of  the 
frequencies  f  =  2,  10,  20  and  30  MHz.  The  results  for  L  =  20  km  are 


- 
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given  in  figure  2  in  comparison  with  the  residue  series  expansion 
result.  It  is  evident  th2t  at  2  MHz  a  step  size  in  excess  of  1  km  is 
adequate  for  convergence,  while  at  30  MHz  the  required  step  size  has 
shrunk  dramatically  to  only  10  m. 

These  results  are  more  conveniently  shown  in  figure  3  which  reprerents 
two  orthogonal  planes  through  the  3-dimensional  space  of  step  size, 
frequency  and  path  length,  showing  the  boundary  of  the  volume  within 
which  the  Dasic  transmission  loss  has  converged  to  within  1  dE. 
Figure  3(a)  indicates  an  approximately  inverse  relationship  oetween  path 
length  and  step  size  for  a  given  frequency,  while  figure  3(b)  indicates 
an  approximately  inverse  square  relationship  between  frequency  and  step 
size.  We  can  therefore  state  an  empirical  convergence  criterion  in  the 
form 


where  X  is  the  wavelength  in  metres,  Ax  is  in  metres  and  L  is  in  km,  and 
the  constant  k  is  yet  to  be  found. 

3.1.2  Effect  of  terrain  gradient 

Real  terrains  exhibit  slopes  which  vary  from  approximately  zero  for  flat 
ground  to  very  large  values  for  almost  sheer  cliff  faces.  Large 
gradients  wori.'  not  normally  be  associated  with  Gaussian  hills  since 
mechanical  stability  considerations  are  likely  to  ensure  that  over 
geological  time-scales  the  slopes  of  hills  would  settle  towards 
equilibrium  values  of  fairly  modest  size,  say  0.5  or  less.  In  spite  of 
this,  we  have  found  it  convenient  to  examine  the  influence  of  terrain 
gradient  over  a  large  range  of  values  by  making  use  of  idealised 
Gaussian  hill  geometries  and  assuming  that  similar  slopes  will  have 
similar  effects  on  the  performance  of  WAGNER,  regardless  of  the  type  of 
topography  involved.  This  approach  follows  the  lead  of  Ott (ref . 2) . 

The  geometry  employed  for  this  study  was  based  on  that  used  for  figure  1 
but  with  the  Gaussian  hill  height  varied  up  to  a  maximum  of  3000  m, 
corresponding  to  a  maximum  gradient  |y'|  ~3,  Figure  4  relates  values 
|y' | ,  and  X  which  just  produce  convergence  with  a  step  size  Ax  =  100  m. 
A  quadratic  curve  which  conveniently  fits  the  data  points  is 
described  by 


X  =  18  (8|y *  I  2  -  2 |y ' |  +  1) 


3.1.3  Effect  i.f  ground  electrical  constants 

It  was  found  that  the  parameters  ground  conductivity  o  and  dielectric 
constant  K  have  a  negligible  effect  on  convergence. 

3.1.4  Convergence  criterion 

We  may  use  equation  (1)  to  write  the  condition  for  convergence  as 
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where  the  function  g  is  identically  unity  (Section  3.1,3)  and  the 
function  f  equals  unity  for  a  slope  of  zero.  For  non-zero  slopes, 
comparison  with  equation  (2)  shows  that  k  =  4.6  and  f  =  8 | y ' |l-2|y' |+1. 
If  we  conservatively  adjust  k  downward  to  k  =  2,  providing  a  safety 
margin  of  about  a  factor  of  2,  we  arrive  at  our  convergence  criterion 
given  by 


Ax  <  2X2L~ 1 (8 1 y ’ | 2  -  2|y’|  +  l)"1 


(3) 
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required  integration  step  size  in  metres 

wavelength  (10  ^  X  $  30C  ra) 

path  length  (0  <  L  $  100  km) 

maximum  value  of  terrain  slope  along  the  path 
(0  <  |y'i  <  3) 


Being  an  empirical  result,  equation  (3)  holds  only  approximately  within 
the  tested  parameter  volume  and  cannot  be  guaranteed  to  hold  elsewhere .- 


3.2  Accuracy 

The  existence  of  mathematical  convergence  is  no  assurance  of  physical 
accuracy.  Therefore  we  varied  parameters  within  the  constraints  of 
equation  (3)  in  order  to  investigate  how  closely  the  convergent  results  of 
WAGNER  agreed  with  basic  transmission  losses  calculated  by  alternative 
methods,  namely  the  residue  series  expansion  for  smooth  Earth  propagation, 
and  the  graphical  knife-edge  technique  for  single-hill  terrains. 

3.2.1  Effect  of  path  length 

It  was  found  that  for  smooth  spherical  Earth  propagation,  WAGNER  always 
ccnverged  to  within  0.1  dB  of  the  value  of  basic  transmission  loss 
predicted  by  the  residue  series  expansion.  This  result  appeared  to  hold 
everywhere  within  the  same  volume  of  parameter  space  as  was  used  for  the 
convergence  study.  In  particular,  path  length  variations  appeared  not 
to  affect  the  accuracy  of  the  final  result,  provided  convergence  had 
been  achieved. 


3.2.2  Effect  of  terrain  gradient  and  frequency 

Figure  5  shows  how  the  accuracy  of  WAGNER  varies  with  frequency  in  the 
case  of  a  £00  m  high  hill  located  5  km  from  the  Tx  antenna  along  a  total 

path  of  15  km,  with  o  =  0.01  S  ra  1  and  K  =  10.  The  convergent  values  ct 
basic  transmission  loss  are  plotted  and  compared  with  the  loss  for 
smooth  spherical  Earth  terrain  and  with  the  addition  of  the  shadow  loss 
predicted  by  the  knife-edge  graphical  technique.  It  is  evident  that 
there  is  agreement  up  to  a  maximum  frequency  of  only  about  2  MHz  and 
chat  above  5  MKz  WAGNER  incorrectly  predicts  a  shadow  gain  with  respect 
to  smooth  spherical  Earth  propagation,  rather  than  a  shadow  loss. 

This  result  and  those  of  similar  calculations  for  hills  up  to  3000  a  in 
height,  corresponding  to  terrain  slopes  |y' |  up  to  about  3.0,  indicate 
that  accuracy  is  governed  by  a  condition  of  the  approximate  form 


ly'i  f  *  i 
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where  f  is  the  frequency  in  MHz.  This  contrasts  strongly  with  the 
"heuristic  uncertainty  principle"  proposed  by  Ott  (see  section  1  above). 
Whereas  Ott's  result  for  |y'|  =  0.2S  predicts  accuracy  up  to  a  maximum 
limiting  frequency  of  40  MHz,  our  equation  (4)  limits  this  maximum 
usable  frequency  to  only  4  MHz. 

3.2.3  Effect  of  ground  electrical  constants 

Figures  6  and  7  present  the  results  of  a  study  of  the  effect  of  ground 
conductivity  a  and  dielectric  constant  K  on  the  accuracy  of  the  integral 
equation  solution.  Figure  6  is  drawn  in  terms  of  shadow  loss  (or  gain) 
versus  frequency  for  different  choices  of  the  ground  constants  and  for  a 
terrain  slope  of  0.25.  It  shows  that  WAGNER  agrees  best  wich  the 
knife-edge  prediction  over  the  1  to  30  MHz  frequency  interval  when  the 
ground  constant  values  are  highest.  In  particular,  for  the  low  values 
of  ground  constants  applicable  to  the  Alice  Springs  region 

(c  =  0.003  S  m  l,  K  =  3),  a  terrain  slope  of  0.25  can  be  accurately 
modelled  up  to  a  maximum  frequency  of  only  2  MHz.  Calculations  for 
additional  terrain  slopes  resulted  in  figure  7,  showing  the  maximum 
frequency  which  can  be  accurately  used  as  a  function  of  the  ground 
constants. 

In  general,  a  and  K  vary  roughly  in  proportion  to  each  other  according 
to  ground  conditions  and  we  find  it  convenient  to  treat  the  product  oK 
as  a  single  key  parameter.  Plotting  its  square  root  as  in  figure  7 
permits  the  data  to  be  represented  on  a  more  compressed  horizontal 
scale.  The  results  of  figure  7  indicate  that  the  maximum  frequency 
which  can  be  used  in  order  to  obtain  accurate  results  is  inversely 
proportional  to  the  terrain  slope  (in  agreement  with  equation  (4))  and 
approximately  proportional  to  the  fourth  root  of  the  product  oK. 

3.2.4  Accuracy  criterion 

The  above  results  lead  to  an  empirical  accuracy  criterion  given  by 

ly’l  f  (oK)"*  <  1.8  (5) 


In  the  Alice  Springs  region  we  may  have  typically  |y'|~0.25, 

o~0.003  S  m  1  and  K~3,  implying  that  accurate  basic  transmission  loss 
results  can  be  obtained  up  to  a  maximum  frequency  of  only  about  2  MHz. 
In  this  case  therefore,  WAGNER  cannot  be  used  to  calculate  basic 
transmission  losses  at  frequencies  in  the  HF  band. 

Alternatively,  equation  (5)  shows  that  at  a  frequency  of  10  MHz  in  the 
Alice  Springs  region,  accuracy  can  be  assured  only  for  terrain  slopes 
less  than  about  0.Q6  or  an  elevation  change  of  60  m  over  a  distance  of 
1  km.  Thus  our  version  of  WAGNER  appears  to  be  limited  to  low 
frequencies  and  gentle  terrains,  if  basic  transmission  loss  is  to  be 
calculated  with  acceptable  accuracy. 

3.3  "Pathological"  behaviour 

While  the  results  described  above  are  disappointing  in  terms  of  possible 
applications  of  WAGNER  in  the  HF  band,  they  nevertheless  indicate  a  smooth 
and  regular  variation  wit.h  respect  to  the  key  parameters.-  However 
occasional  sharp  departures  from  this  predictable  behaviour  have  been 
observed.  For  example,  in  the  case  of  smooth  spherical  earth  propagation 

at  10  MHz,  with  a  -  0.003  S  m  *,  K  =  3,  L  =  100  km,  and  Ax  =  100  m,  the 
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calculated  basic  transmission  loss  appeared  to  be  converging  as  expected. 
However  further  reduction  of  the  step  size  to  50  m  surprisingly  produced  a 
sharply  divergent  result  for  distances  in  excess  of  50  km.  Additirnal 
reduction  in  step  si2e  then  led  to  a  recovery  with  convergence  finally 
being  achieved  for  Ax  =  20  ra. 

The  user  is  warned  of  the  possibility  of  occasional  rare  occurrences  of 
this  type.  One  safeguard  is  to  always  examine  the  trend  of  the  basic 
transmission  loss  as  a  function  of  distance  along  a  path,  looking  for 
irregular  behaviour. 


4.  DISCUSSION 

The  results  of  Section  3  above  show  that  even  when  a  calculation  by  WAGNER  of 
basic  transmission  loss  has  converged  to  a  stable  value,  the  result  remains 
subject  to  a  severe  physical  accuracy  constraint  which  renders  our  version  of 
WAGNER  almost  useless  in  the  KF  band.  This  is  particularly  so  when  the  ground 
electrical  constants  have  low  values,  as  in  the  case  of  the  arid  desert 
environment  surrounding  Alice  Springs. 

However  in  many  applications  the  total  basic  transmission  loss  may  not  be  the 
required  variable.  It  may  suffice  to  calculate  only  the  shadow  loss  (relative 
to  smooth  spherical  Earth  propagation)  due  to  one  or  more  terrain 
irregularities..  Since  shadow  loss  depends  mainly  on  topography  and  very 
little  on  the  electrical  constants  of  the  terrain,  it  may  be  possible  to 
calculate  shadow  loss  by  using  "dummy"  values  of  conductivity  and  dielectric 
constant  which  are  chosen  to  be  high  enough  to  ensure  numerical  accuracy 
(equation  (5)).  The  program  could  be  run  first  for  the  terrain  under 
consideration  and  then  for  a  smooth  spherical  Earth,  and  the  first  result 
divided  by  the  second  in  order  to  obtain  the  shadow  loss..  If  shadow  loss 
dcpar.is  only  on  the  frequency  and  on  the  geometrical  dimensions  of  the 
obstacle,  a  reliable  result  should  be  obtained,  with  the  dummy  ground 
constants  being  self-cancelling.  For  the  case  of  a  single  obstacle  this 
result  should  agree  well  with  the  graphical  estimate  using  reference  8. 
In  the  case  of  arbitrarily  complex  terrains  WAGNER  alone  possesses  the 
potential  to  produce  believable  shadow  loss  results,  subject  to  the 
constraining  equation  (5). 

We  are  able  to  demonstrate  the  usefulness  of  this  approach.  Figure  8  compares 
the  shadow  loss  predicted  by  the  knife-edge  technique  with  two  results  from 
WAGNER,  firstly  a  single  Gaussian  hill  approximation  to  the  terrain  and 
secondly  the  actual  terrain  interpolated  into  intervals  equal  in  size  to  the 
specified  integration  step  size  Ax..  The  WAGNER  calculations  were  carried  out 

for  a  dummy  conductivity  of  50  S  m  '  and  a  dummy  dielectric  constant  of  500. 
Also  shown  are  actual  shadow  loss  data  measured  over  this  path.  The 
calculations  using  the  actual  terrain  profile  agree  best  with  the  observations 
at  the  lower  frequencies  but  tend  to  overestimate  the  loss  above  15  MHz.  The 
Gaussian  profile  predicts  a  frequency  variation  of  shadow  loss  which  is  not 
matched  by  the  rather  flat  variation  of  the  observations.  The  graphical 
knife-edge  estimates  also  fall  away  with  frequency  rather  more  rapidly  than  is 
actually  measured,  but  are  nevertheless  in  fair  agreement.  Application  of  the 
technique  to  other  paths  has  also  been  generally  successful(ref .8) . 

The  above  rusults  indicate  that  shadow  loss  in  the  HF  band  may  be  reliably 
obtained  from  WAGNF.R  provided  firstly  that  the  convergence  criterion  is 
satisfied,  and  secondly  that  the  accuracy  criterion  is  satisfied  by  using  high 
dummy  values  of  ground  conductivity  and  dielectric  constant.  That  the  method 
proves  to  be  satisfactory  relies  on  shadow  loss  being  largely  independent  of 
ground  electrical  constants. 
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5 .  CONCLUSION 


In  the  absence  of  any  reporced  theoretical  analysis,  we  have  empirically 
derived  heuristic  criteria  which  govern  the  convergence  behaviour  and  physical 
accuracy  of  Ott's  integral  equation  approach  to  the  calculation  of  ground  wave 
propagation  losses  over  irregular  and  inhomogeneous  terrain.  We  have  found 
that  for  irregular  terrain  the  existence  of  mathematical  convergence  does  not 
necessarily  imply  a  physically  accurate  calculation. 


The  computer  program  WAGNER  which  implements  the  integral  equation  calculation 
has  been  tested  for  homogeneous  ground  by  varying  selected  key  parameters  over 
intervals  pertinent  to  ground  wive  propagation  work  being  carried  out  in 
conjunction  with  HF  radar  studies.  These  parameters  and  their  test  intervals 
are  as  follows: 


Integration  step  size  (uniform)  Ax 
Frequency  f 

Path  length  L 

Terrain  maximum  slope  iy'l 

Product  of  ground  conductivity  (o  K) 
and  dielectric  constant 


10  to  1000  m 
1  to  30  MHz 
10  to  100  km 
0  to  3.0 

10'J  to  10‘  S  m" 1 


The  following  heuristic  criteria  were  found  to  govern  the  mathematical 
convergence  and  physical  accuracy  of  the  calculated  basic  transmission  less. 


Mathematical  convergence: 


Ax  <  a).2!.- 1  (8 |y *  |  2  -2|y*  |  +  l)_l 


Physical  accuracy: 


|y’|  f  (oK)"*  <  1.8 


It  is  stressed  that  the  empirical  origin  of  these  expressions  means  that  they 
are  not  unique  and  their  validity  cannot  be  guaranteed  outside  the  tested 
parameter  ranges  specified  above. 

While  our  results  confirm  the  usefulness  of  WAGNER  m  calculating  basic 
transmission  loss  near  a  frequency  of  1  MHz,  they  demonstrate  serious  accuracy 
difficulties  at  higher  frequencies.  In  the  case  f  =  10  MHz,  |y'|  =  0.25, 

a  =  0.003  S  m  l,  K  =  3  and  L  =  20  km,  mathematical  convergence  may  be  achieved 
protided  the  integration  step  size  Arc  is  reduced  to  about  100  m.  On  the 
HF  Radar  Division  VAX  8200  running  under  the  operating  system  VMS  version  4.4 
and  using  VAX  FORTRAN  version  4.5,  the  execution  time  for  this  calculation  is 
approximately  three  minutes.  Convergence  is  therefore  not  a  difficulty  except 
possibly  at  the  high  extremes  of  terrain  slope  and  patn  length  when  small  step 
sizes  and  hence  a  large  number  of  steps  are  required.  Since  execution  time 
varies  as  the  square  of  the  number  of  integration  steps  such  a  calculation  rosy 
run  for  many  hours,  but  convergence  will  ultimately  be  reached.  On  the  other 
hand,  physical  accuracy  appears  as  a  much  more  serious  problem  and  basic 
transmission  loss  over  the  above  path  can  be  modelled  by  WAGNER  to  an  accuracy 
of  better  than  3  dB  only  for  frequencies  below  2  MHz.  The  presence  of  hills 
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of  only  modest  size  therefore  effectively  precludes  the  use  of  WAGNER  for  the 
calculation  of  basic  transmission  loss  in  the  HF  band.  This  result  does  not 
appear  to  have  been  spelt  out  by  the  authors  of  previous  reports  dealing  with 
WAGNER. 

If  it  is  desired  only  to  determine  the  relative  shadow  loss  at  HF  due  to 
irregular  terrain  along  the  path,  it  is  feasible  to  perform  the  integral 
equation  calculation  using  unrealistically  high  ground  electrical  constants  in 
order  to  guarantee  physical  accuracy  and  then  compare  the  result  with  the 
smooth  spherical  Earth  calculation  carried  out  with  the  same  parameters.  The 
dummy  electrical  constants  effectively  cancel  themselves  out  in  the  process 
and  therefore  do  not  affect  the  derivation  of  the  shadow  loss.  The 
demonstrated  viability  of  this  technique  permits  WAGNER  to  be  used  at  HF  for 
arbitrary  multiple  hills,  which  are  not  easily  treated  by  the  graphical 
knife-edge  technique. 

Importantly,  our  condition  on  physical  accuracy  is  in  disagreement  with  the 
condition  given  by  Ott(ref.l).  Whereas  Ott's  condition  implies  that  terrain 
slopes  of  the  order  of  1  in  4  can  be  treated  throughout  the  HF  band  up  to  a 
maximum  frequency  of  40  MHz,  we  find  that  our  upper  frequency  limit  is 

actually  4  MHz,  assuming  a  =  0.01  S  m  1  and  K  =  10;  the  HF  band  is  therefore 
effectively  excluded  from  calculations  of  basic  transmission  loss  unless 
larger  values  of  a  and  K  can  be  employed. 

Finally,  we  emphasize  that  the  results  presented  in  this  report  were  obtained 
using  our  double-precision  version  of  WAGNER,  locally  known  as  W5D.  It  is 
possible  that  more  robust  versions  of  WAGNER  have  evolved  elsewhere.  Every 
effort  has  been  made  to  acquire  versions  currently  in  use  in  other  parts  of 
the  world  but  without  success.  However  this  effort  is  continuing.  The 
authors  will  be  happy  to  supply  a  copy  of  W5D  upon  request. 
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FORTRAN  program  'WAGNER' .  The  potential  vaiue  of  this  integral  equation 
technique  is  its  purported  ability  to  model  propagation  losses  over 
'  two-dimensional  paths  through  Irregular  and  inhomogeneous  terrain.  However 
the  numerical  Integration  algorithm  can  be  affected  by  poor  convergence  and 
inaccuracy,  particularly  at  frequencies  above  the  MF  band.  We  have 
investigated  these  properties  in  terms  of  selected  key  parameters  and  have 
derived  empirical  criteria  governing  convergence  and  physical  accuracy. 

We  find  that  calculations  of  basic  transmission  loss  using  WAGNER  are 
subject  to  significant  limitations  at  HF  frequencies  as  the  terrain  becomes 
irregular  and  as  the  conductivity  and  dielectric  constant  of  the  ground 
both  become  smaller.  On  the  other  hand,  the  shadow  loss  due  to  an  obstacle 
may  be  reliably  calculated  by  WAGNER  on  the  assumption  that  shadow  loss  is 
independent  of  ground  electrical  constants,  thus  permitting  the  use  of 
|  dummy  values  of  conductivity  and  dielectric  constant  to  ensure  that  the 
|  physical  accuracy  criterion  is  satisfied.  / 


